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Abstract 

We consider the 2D Ising lattice gas in a square of side L with free 
boundary conditions, temperature below the critical one and particle 
density slightly above the density of the vapor phase. Typical con- 
figurations consist of a quarter of a Wulff droplet of the liquid phase 
centered at one of the corners of the given square. We then intro- 
duced a reversible Markovian spin exchange dynamics, also known as 
Kawasaki dynamics, on the configuration space and we discuss the 
heuristics of the transition of a bubble of the liquid phase from one 
corner to another. We then present some numerical evidence sug- 
gesting that the preferred mechanism to make the transition is via 
evaporation of the original bubble and simultaneous reconstruction of 
a new bubble around a new corner. 
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1 Introduction 



The problem of computing the relaxation time of stochastic Monte Carlo al- 
gorithms for models of lattice classical spin models has attracted in the last 
years considerable attention and many new rigorous techniques have been 
developed giving rise to nice progresses in probability theory and statistical 
mechanics. If for simplicity we confine ourselves to ±1 (or 0, 1 in the lattice 
gas picture) spins, the two most studied random dynamics have been non- 
conservative Glauber type algorithms, in which a spin at a time flips its value 
with a rate satisfying the detailed balance condition w.r.t. the grand canon- 
ical Gibbs measure, and conservative Kawasaki dynamics in which nearest 
neighbors spins exchange their values with a rate satisfying the detailed bal- 
ance condition w.r.t. the canonical Gibbs measure. 

A warning sign here is in order. Most of the rigorous results presented 
below have been obtained in the "continuous time" setting, namely when the 
stochastic dynamics is a continuous time Markov chain and the time units 
are such that, in a unitary time interval on average every spin tries to change 
its value once. In the physics literature, but also in the Markov chain Monte 
Carlo approach to computational problems and image reconstruction, the 
discrete time setting in which at each time step only one dynamical variable 
is updated is more familiar. In most cases the simple rule to translate results 
from one case to another is to multiply by the appropriate change of scale, 
but sometimes subtle problems may appear in the discrete setting (see e.g. 

i)- 

Let us go back to a quick review of the main results. For Glauber dy- 
namics in the one phase region the general picture is relatively clear for a 
wide class of models and the conclusion is that equilibrium is reached expo- 
nentially fast provided certain mixing assumption are satisfied by the grand 
canonical Gibbs measure (see e.g. and references therein). In the con- 
servative case, instead, the basic results [||, and more recently 0, state 
that, at high temperature, the relaxation time in a box of side L grows like 

(diffusive scaling). 

A natural question arises as to what happens when the thermodynamic 
parameters (e.g inverse temperature and the external magnetic field in the 
Glauber case or inverse temperature and particle density in the conservative 
case) are such that we do have a phase transition in the thermodynamic 
limit. 

Let us start with the Glauber case and, to be concrete, let us consider 
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the usual Ising model in a two dimensional box V of side L, without external 
field and inverse temperature (3 larger than the critical value (3c- With free 
boundary conditions the picture of the relaxation behavior to the Gibbs 
equilibrium measure that comes out is the following. The system first relaxes 
rather rapidly to one of the two phases and then it creates, via a large 
fiuctuation, a thin layer of the opposite phase along one of the sides of V. 
Such a process requires an average time of the order of exp(/9r^L) where 
denotes the surface tension in the direction of one of the coordinate axes. 
After that, the opposite phase invades the whole system by moving, in a 
much shorter time scale, the interface to the side opposite to the initial one 
and equilibrium is finally reached. The time required for this final process 
can be computed to be of the order of at least in the SOS approximation 
(see i). 

Let us now turn to the much more involved conservative case in the same 
setting (low temperature and free boundary conditions). Here results are 
much less precise (see 0). 

Let us denote by N the total number of particles in V and let us assume 
that p e (pl,p+), where p, denote the actual density and the density 
of the liquid/vapor phases respectively. With these assumptions it can be 
shown that the relaxation time is again exponentially large in L but, con- 
trary to the conservative case, almost nothing is known rigorously about the 
exact constant in the exponential and about the physical mechanisms behind 
equilibration. 

That the relaxation time is very large can be easily understood (but 
painfully proved) by the following reasoning. Assume for simplicity that p 
is slightly above p*_ in such a way that, typically, the gas at equilibrium 
shows phase segregation between a (macroscopically) small bubble of liquid 
immersed in vapor. Because of the free boundary conditions the bubble 
prefers to sit around one of the corners of V. Clearly, in order to reach 
equilibrium, the system started as above has to move the bubble to the other 
corners but, in doing that, it is forced to explore a region of the phase space of 
exponentially (in L) small canonical probability (and that is the hard part in 
a rigorous approach)). In other words it seems clear that the slowest mode in 
the dynamics comes from the corner to corner droplet transport. A precise 
understanding of the above motion appears therefore to be an interesting 
problem on its own and almost unavoidable if one wants to quantify precisely 
the relaxation time. A first numerical attempt in this direction represents 
the main goal of this note, which is in turn based on the work contained in 
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i- 

We conclude by observing that, if boundary conditions are changed to e.g 
plus or periodic, then the whole scenario changes drastically and one may 
argue that the slowest mode of the system is associated to the random walk 
motion of the center of gravity of the unique Wulff bubble of the liquid phase. 
If that is true then simple heuristics shows that the relaxation time should 
now grow as Li^ (in continuous time units) or in discrete time steps. 

2 Model and Notation 

We define here our model and notation. 

2.1 Equilibrium Probability Measures 

We consider the two dimensional lattice 7? with sites x = (xi,X2). By Ql 
we denote the square of all x = {xi, X2) G such that Xj G {0, . . . , L — 1}, 
i = 1,2. The edges of are those e = [x,y] with x,y nearest neighbors 
in Z^. Given C Z^, we denote by Sy the set of all edges such that both 
endpoints are in V. 

The configuration space. Given C Z^, our configuration space is fly ~ 
S"^, where S = { — 1, 1}. Sometimes the lattice gas point of view will be more 
convenient, so we also consider the space Qly = {0, 1}^ and its natural one- 
to-one correspondence with Q. We define the unnormalized magnetization 
My '. ^ Z and the number of particles Ny '■ ^'y ^— N as 

My{a) = Y.a{x), Nyiv) = J2r^ix) (1) 

while the normalized magnetization is given by my = My /\V\. 

The interaction and the Gibbs measures. Given a finite set C Z^ we 
define the Hamiltonian with free boundary condition Hy : Qy 1— R by 

Hyia) = - 5: aix)aiy). (2) 

The corresponding (finite volume) Gibbs measure is given by 

f,^y{a) = {Z^)-'exp[-(3Hy{a)] , (3) 
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where Zy is the proper normahzation factor called partition function. We 

denote by m*{j3) and j3f. the spontaneous magnetization and inverse critical 
temperature of the two dimensional Ising model respectively. It is well known 
that (3c = (1/2) log(l + ^/2). The density of the hquid and the vapor phase, 
denoted respectively by and p*_{P), can be expressed as |(1 =bm*). 

We also introduce the canonical Gibbs measure defined as 

4,N = l^vi-\Nv = N) N e{Q,l,...,\V\} , (4) 
where Ny is the number of particles in V. 



3 The Kawasaki Dynamics 

Here we define the relevant Markovian dynamics, reversible with respect to 
the canonical Gibbs measure, that will be used and analyzed in the rest of 
this paper. 

Given V G Z'^, b ^ Sy, and a particle configuration ri G Q'y (equivalent 
to a spin configuration a G 0,v), let be the configuration obtained from rj 
by interchanging the values of the rj variables at the end points of the bond 
b. The energy difference between the two configurations is given by 

A,Hviv) = Hviv") - Hviv) ■ (5) 

II ^ T] we call b a broken bond and we say it belongs to B^. With the 
above notation the Kawasaki dynamics with Metropolis transition probability 
matrix is given by: 

f ^e-^°^^(0'^i''^^('')) if 77' = rf and 6 G , 
W{r^, V) = 1 - W{r^, rj") if 77' = 77 , (6) 

. otherwise . 

Clearly, for each G [1, |V^| — 1], 14^ describes an ergodic Markov chain on 
the configuration space fl'vN which consists of all particle configurations with 
particle number equal to N. In particular, since W satisfies the detailed 
balance condition with respect to the canonical measure i^yj^ latter 
coincides with the unique invariant measure of the chain. 
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4 Corner to Corner Matter Transport 



In this section we will first try to define our field of investigation by start- 
ing with a discussion of the typical configurations for the canonical measure 
of the Ising model below the critical temperature with free boundary con- 
ditions. We will then analyze some possible mechanisms for the corner to 
corner bubble transition, and the typical time scales that govern the process. 
Subsequently we will discuss the details of our numerical simulations and 
give some hints about the updating algorithm. Finally we will present the 
main results of this note. 



4.1 Heuristics 

Let V = Ql, let p G {p*_,p\), where have been defined above, and let 
N = \_pL'^\ be the total number of particles. 

For the above situation, it is useful to recall first the shape of the typical 
configurations of the canonical Ising Gibbs measure with free b.c. when the 
temperature is below the critical value and L is very large. 

Let rUp = 2p — 1 be the usual magnetization associated to the given 



particle density. Then, as discussed in [|Tol and 0, there exists < mi{(3) < 
m*{l3) such that 

1. (i) if rup G {—mi, mi) then the typical configurations show phase segre- 
gation between a high density (~ p+) region and a low density (~ pi) 
region that are roughly two horizontal (or vertical) rectangles of appro- 
priate area separated by an horizontal (or vertical) interface of length 

2. (ii) if mp G (-m*(/?), -mi{(3)]VJ[mi{(3),m*{(3)) then the typical config- 
urations show phase segregation between a high density (~ p\) region 
and a low density (~ pi) region, the smaller of which is a quarter of 
a Wulff shape (see [|12|) of appropriate area and it is centered in one 
of the four vertices of Ql. For the reader convenience we recall that a 
Wulff shape is, at very low T, very similar to a perfect square. 

What is important for us is that in both cases the typical configurations 
of the canonical measure show a discrete symmetry described by rotations of 
= 0, 1 . . . around the center of A. As a consequence, if the dynamics 
starts from one typical configuration for which for example the particles form 
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a cluster in a corner (this is one of the situations we have described before), 
then, in order to reach equihbrium, it must necessarily cross an unlikely 
region in the configuration space. One can show that the canonical proba- 
bility of such region is exponentially small in L . Therefore a bottleneck is 
present in the configuration space, and the relaxation time is exponentially 
large in L 

In the sequel of this note we will always assume to be in the second 
between the two scenario's described above, and in particular we will assume 
that m and p are respectively larger but quite close to the values — m*(/3) 
and p*_{P)- In other words we are considering a situation in which typically 
the particles form a small (on a macroscopic scale) cluster of a specific shape 
around one of the four corners of A. Clearly we do not mean that all the 
particles that are in A are in the bubble, and we do not mean that the 
bubble does not have empty holes in its interior. What we mean exactly is 
that, with high probability, there exists a macroscopic region with a precise 
shape around one of the corners where the particle density is very close to the 
density of the liquid phase, p\ (corresponding to the Onsager magnetization 
density) . 

Next we analyze some possible mechanisms for the corner to corner bubble 
transition, namely the process that moves the liquid bubble from one corner 
to a different one. 

We can see two main mechanisms that could intervene. 

The first one involves a deformation of the initial Wulff droplet, that for 
simplicity we will assume here and in the following to be a square of side B 
(where pL^ = p\B^ + p*_{L'^ — B"^)) and with one vertex sitting on the left 
lower corner of A, into a rectangle. We call this case sliding. We can describe 
it as follows: we assume that the particles always (that is also during the 
transition to one of the other corners) form one compact cluster (apart from 
the usual small fluctuations) that, based on energetic considerations, we can 
assume to be a generic rectangle IZ of sides Bi and B2. Energetically, because 
of the free boundary conditions, at equilibrium it will find convenient to be 
a square of side B attached to two sides of A. Under the Kawasaki dynamics 
we can imagine that such square gets deformed, by some sort of matter 
transport along the boundary, until it reaches the shape of a rectangle TZ* of 
sides Bi > B2 [Bi being the horizontal side) with the left and bottom sides 
attached to the boundary of A. At this point the rectangle TZ* begins to slide 
along e.g. the bottom side of A until it reaches the opposite corner where it 
deforms back to the original square (see figure (|l])). The energy barrier AH 
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Figure 1: The sliding motion of the bubble, from the left to the right bottom 
corner. 



one has to cross in this process is clearly of the order of 2{Bi + 2B2) — 2(25) 
at least at very low temperature. If we now optimize over Bi and B2 under 
the constraint that Bi x B2 ^ B^ we get Bi = y/2B and B2 = ;^-B. Thus 

AH = 4(\/2 — l)B and therefore wc expect that the average time to see a 
sliding transition to be of the order t sliding ~ exp(/5AiJ). 

The second mechanism that one can imagine is what wc call evaporation. 
Particles individually separate from the liquid bubble (they evaporate) and 
start to perform a sort of random walk in the vapor phase (typically along the 
boundaries, by energetic reasons). Clearly the excess particles in the vapor 
phase prefer, as soon as they can, to gather together around one of nearby 
corners, and grow there a new liquid bubble. Typically these trial bubbles do 
not achieve a macroscopic volume, but die much before reaching this stage 
and their particles end up to go back to their mother bubble. Only rare fluc- 
tuations lead to more than half of the particles clustering together: in these 
cases typically the new bubble will form in the selected corner, and evap- 
oration will continue until all the original liquid bubble has been reformed 
around the new corner (clearly particles of the bubble get interchanged at 
all moments with particles that are in equilibrium in the vapor phase). It 
is not difficult to check that, at least to leading order, the energy barrier is 
crossed when the two groups of particles in the two different corners have 
the same volume (and same WulfF shape). At very low temperature we get 
immediately that AiJ — 4(\/2 — 1)B, exactly as in the former case. 

We conclude this part by stressing that all the above reasonings were 
based only on energy barrier considerations and that we have never taken 
into account cntropic contributions. The latter may play an important role 
in selecting one mechanism instead of another. Moreover prefactors in the 
typical time scales of each process may also be relevant and in that case a 
more detailed analysis would be required. 
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4.2 Numerical Simulations 



Our experimental setting can be described as follows: we use a square two 
dimensional lattice with free boundary conditions. We have investigated 
lattice sizes going from L = 20 to L = 30. The inverse temperature j3 has 
been assigned the two values (3 = 0.7 and (3 = 1.05 in different cases (this is 
of the order of twice the value of the inverse critical temperature /3c). Finally 
we have chosen the number of particles in the range 25 to 36 (corresponding 
to initial conditions with a bubble of 5 • 5 or 6 ■ 6 particles in a corner of the 
square lattice). Notice that in our range of temperatures the density of the 
vapor p*_ lies between 5 10~^ and 2 10"'^ for the two different temperature 
values: in other words we are in extreme conditions. 

The updating algorithm is of course based on the transition matrix (^, 
and goes as follows. We consider the particle configuration at time t, r]{t). We 
select a broken bond (that is a bond in the set S^(t) ) at random, with uniform 
probability. We compute the energy difference ^}yHy{rj{t)) defined in (^: if it 
is negative we accept the update proposal, and set rj{t+l) = rfit). Otherwise 
we accept the update proposal with probability exp{—[3^bHv{r]{t))) , and 
refuse it otherwise, setting in this case rj{t + 1) = rj{t). 

This procedure, while appealing from the point of view of the computa- 
tional efficiency, and satisfactory from the physical point of view, does not 
satisfy detailed balance, since the number of broken bonds is not a conserved 
quantity (see for instance the discussion in [llj). In our working conditions, 
however, this effect is very small and as far as many issues are concerned, 
irrelevant. In fact, as we have already explained, we will mostly focus on the 
mechanisms on which the bubble transition is based: since the small viola- 
tion of detailed balance amounts to watching a movie that runs at slightly 
variable speed, spatial phenomena (like for example which is the transition 
path of the bubble or which is its typical spread during the transition) are 
described in an exact manner. Scaling arguments (for example the scaling of 
the transition time with (3) could as matter of principle be sensitive to the 
violation of detailed balance: we believe however that only the prefactors 
will be affected, and all of the many tests of universality that we have done 
on our simulations do confirm this point of view. 

We still have to discuss the criterium by which we define a bubble tran- 
sition. We have indeed used two different criteria and checked that in the 
two cases we get consistent results. In the first scheme we define four boxes, 
with a vertex in one of the four corners, with an area equal to the one of 
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the original bubble (that starts, let us say, from the top left corner). We 
then define a transition as the event where the center of mass of the particles 
that constitute the bubble enters a new box. In the second scheme we define 
a transition as the event where the | of the particles have left the initial 
box: since we expect to have a transition when half of the particles have left, 
we are confident that 75 percent is a safe signal for a transition. Numerical 
simulations will confirm the coincidence of the two criteria. 

5 Results and Discussion 

In figure (0) we show the number of particles (modulo a corner dependent 
offset needed for making the figure readable) at time t near the 4 different 
corners of the lattice during the course of the dynamics. More precisely we 
have defined square boxes of size with a vertex in each of the four corners, 
where is the total number of particles (25 in the present case) and for 
each box we have computed the time history of the number of particles inside 
the box, under the condition that at time t = all the particles are in the 
leftmost bottom box. The value for the lowest curve (representing one given 
corner) is exactly the value of the number of particles in the corner, while the 
other 3 curves have an offset of, respectively, 30, 60 and 90 for the 3 different 
corners. The transitions are always abrupt, and the particles spend the quasi 
totality of their time confined in one corner. As one expects on theoretical 
grounds the typical time scale on which a transition occurs is much shorter 
than the time one has to wait to see the transition. The system is waiting for 
a critical fluctuation making a lot of unsuccessful attempts (the small spikes 
towards the bottom of the figure). 

In figure (^ instead we have computed the "variance" of the center of 
mass of the particles. We define 

^'W = iE(^^W-^^'^™nt))' , (7) 

i 

where Xi{t) cirG, clS usual, the positions of the particles on the two dimensional 
square lattice, and x^'^"^\t) is the center of mass of the particles at time t. 
We compare it to the two extreme situations, where the particles form a 
compact blob (the lowest straight line), and where the particles are divided 
in two, equal sized, compact blobs in two adjacent corners (see the previous 
discussion in §3.1). The time scale of this figure (in arbitrary units again) is 
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Figure 2: The number of particles close to each of the four corners as a 
function of time. See the text for further details. 



much smaller than the one of the previous figure, and basically gives the time 
of a single transition. In other words we have analyzed the time behavior 
of the above variance precisely during the short (in the time scale of figure 
(0)) time interval in which the transition takes place. The outcome is that 
the particles go from a compact blob (in a corner) to a compact blob (in a 
different corner) passing through a situation where they are basically divided 
in two groups of similar size. 

The second part of our analysis considers scaling behaviors of the (aver- 
age) transition times T(/3,L). We look separately at the dependence of T 
over f3 and over the size L. In figure (^ we look at the f3 dependence over 
the range [0.65, 1.1] for a fixed side L = 20 (and 25 particles). Our best fit 
to the form 

T{p) = A e^''^ (8) 

is very good and gives AH ~ 17.8 (and A ~ 0.001). 

This is in very good agreement with the heuristic calculation of the energy 
barrier assuming that the transition occur via the mechanism we have called 
evaporation, that suggests that AH is of order 20. We illustrate the mecha- 
nism in figure (1). We start with 25 particles packed in the left down-most 
corner: here the surface includes 10 broken links. Now in the intermediate 
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Figure 3: Variance of the center of mass as a function of the time. 
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Figure 4: Transition time as a function of f3. 
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Figure 5: The dynamical mechanism we call evaporation. See the text for 
further details. 

situation of figure that is a typical intermediate particle configuration, 
we have 19 broken links. Clearly it would be less expensive for the lonely par- 
ticle to travel along the side of the lattice, but the configuration we show has 
a higher entropy and we observe these kind of processes with high frequency. 
Since AH is equal to twice the difference in the number of broken links, in 
this case AH = 18. This is only an estimate, but gives the correct order 
of magnitude. Notice that AH = 18 is quite far from the infinite volume 
saddle we have discussed before, that would give AH = 4(^2 - 1)B ~ 8: 
this is due to finite size effects (both B and L are finite). We have verified 
that when we change the number of particles (for example from 25 to 36) or 
the definition of a transition from corner to corner we continue to find the 
expected scaling behavior. 

6 Conclusions 

In this paper we have investigated, mainly numerically, the Kawasaki dynam- 
ics ( a particular conservative particle-hole exchange) for the low temperature 
(/? ~ 2/?c) Ising lattice gas in a box with free boundary conditions and side 
L ^ 30. The number of particles N has been chosen in a such a way that 
phase segregation occurs (typically ~ 25). With the above choice of the 
thermodynamic parameters the typical equilibrium configurations consist of 
a (small compared to the total volume) square-like droplet of particles which 



13 



sits in one of the four corners because of the chosen boundary conditions and 
a rarefied gas elsewhere. Under the Kawasaki dynamics the droplet of par- 
ticles make rare and mostly unsuccessful attempts to migrate to one of the 
empty corners. Based on energetic considerations alone, thus neglecting en- 
tropic contributions, we have envisaged two main best possible mechanisms 
for the migration to take place, that we have named sliding (the bubble of 
particles gets deformed into a rectangle, slides along one side of the box until 
it reaches the opposite corner and finally reconstructs the optimal square- 
like shape) and evaporation (the particles in the bubble evaporate into the 
rarefied gas and recollect together to a form a new bubble in another cor- 
ner) respectively. We have then performed intensive numerical investigation 
to check whether sliding or evaporation is the prefered mechanism and the 
scaling behavior with the inverse temperature f3 of the (average) transition 
time. Under various different definitions of the bubble transition our results 
consistently indicate that evaporation is the dominant effect and that the 
scahng law for the average transition time T{P) is of the form 

r(/3) = A e^^"" , (9) 

where the energy barrier 5 H is in very good agreement with simple energetic 
considerations. 
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